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A broad class of dark energy models, which have been proposed in attempts at solving the 
cosmological constant problems, predict a late time variation of the equation of state with redshift. 
The variation occurs as a scalar field picks up speed on its way to negative values of the potential. 
The negative potential energy eventually turns the expansion into contraction and the local universe 
undergoes a big crunch. In this paper we show that cross-correlations of the CMB anisotropy 
and matter distribution, in combination with other cosmological data, can be used to forecast the 
imminence of such cosmic doomsday. 
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I. INTRODUCTION 

Observational cosmology has yielded several surprises, 
of which the most perplexing is the discovery of a smooth 
dark energy (DE) component which has come to domi- 
nate the universe at recent epochs, causing cosmic ac- 
celeration P, 0]. The nature of this component is still 
a matter of speculation, and a very important challenge 
for the coming years will be to determine its origin and 
physical properties. Recently, several groups J, , 4.i ,5 j ^ lII 
have reported positive results for the cross-correlation 
between the CMB power spectrum and that of different 
LSS surveys, providing further evidence for the existence 
of DE. In this paper, we shall try and illustrate how such 
cross-correlation may help in unveiling some of the prop- 
erties of DE, focusing on the observational signatures of 
a model with a time dependent DE equation of state. 

The simplest interpretation of the dark energy is in 
terms of a cosmological constant, with equation of state 
po = —pD- The cosmological constant, however, raises 
two puzzles of its own. First, there is a glaring discrep- 
ancy between the observed dark energy density po and 
the huge values of the cosmological constant suggested 
by particle physics models. Second, the observed is 
comparable to the matter density pm- This is the notori- 
ous time coincidence problem: why do we happen to live 
at the epoch when the dark energy starts dominating? 

It has long been suggested that both puzzles may 
find a natural explanation through anthropic selection 
effects, in scenarios where po is a random variable, 
taking different values in different parts of the uni- 
verse. The proposed selection mechanism is very simple 
H H m El El m [13 . The growth of density fluctua- 
tions leading to galaxy formation effectively stops when 
Pd comes to dominate over the matter density. In re- 
gions where is greater, it will dominate earlier, and 
thus there will be fewer galaxies (and therefore fewer ob- 
servers). A typical observer should then expect to find 
herself in a region where pu dominates at about the 
epoch of galaxy formation (which is close to the present 



time). Much larger values of po would yield no galax- 
ies at all, while much smaller values are unlikely due to 
the smallness of the corresponding range of po , assuming 
that all values of po are a priori equally likely [Tllll5l|. 

A simple implementation of this idea is obtained by 
assuming that the dark energy is due to a scalar field 
with a very flat potential V^(0) J^fi]- The values of 
are randomized by quantum fluctuations during inflation, 
resulting in a variation of </) with a characteristic scale 
much greater than the present Hubble radius. Galaxy 
formation is possible only in regions where V{(j}) is in a 
narrow range near = 0. One expects that the potential 
in this range is well approximated by a linear function 




F(0) = (1) 

where s = —V'{0) and we have set = at y = 0. The 
slope s should be sufficiently small, so that the variation 
of (p is not fast on the present Hubble scale. Quantita- 
tively, this can be expressed as the slow roll condition, 

s < 3H^Mp, (2) 

where Mp = I/vottG is the reduced Planck mass, G is 
Newton's gravitational constant, and Hq is the present 
Hubble expansion rate. 

In models with a single DE field, and in the absence of 
ad-hoc adjustments, it has been argued in E3 that the 
slow roll condition (0) is likely to be satisfied by excess, 
by many orders of magnitude, rather than marginally. 
In this case, remains nearly constant on the Hubble 
scale, and the effective equation of state for the dark 
energy is w = pd/pd ~ ^1, with a very high accuracy. 
However, a different situation may be expected in multi- 
field models, where the slope of the potential s is itself a 
random variable [T^ |2]| (the role of the field (p in multi- 
field models is played by the variable in the direction of 
Vy(0a) in the field space.) The observed value of the 
slope may then be influenced by anthropic selection. A 
very large slope would cause a big crunch much before 
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any observers can develop. Hence, in cases where the 
prior distribution favors large s, the most probable values 
of the slope would be the ones for which the slow roll 
condition ^ is only marginally satisfied (see also 

We thus have some motivation to consider a model 
where the dark energy is due to a scalar field with a lin- 
ear potential and a slope s marginally satisfying the 
slow roll condition |5J). A marginal value of the slope 
implies that the big crunch is imminent in about 10 bil- 
lion years from now. The model can therefore be called 
a "doomsday model" . A salient feature of this scenario 
is that the equation of state of dark energy changes sig- 
nificantly at low redshift, when the correlation between 
the large scale structure evolution and the CMB temper- 
ature anisotropics develop. Hence, we may expect that 
the analysis of such cross correlation may reveal a time 
varying equation of state po — 'w{z)pD, where z is the 
redshift. This will be the subject of the present paper. 

Prospective constraints on cosmic doomsday based on 
future determinations of the dimming of distant super- 
nova were discussed in Ref. 22]. The analysis shows that 
SNIa observations, in combination with CMB and weak 
lensing data, have an impressive potential for constrain- 
ing the equation of state parameter. However, the con- 
straints reported in still show a considerable degen- 
eracy amongst models with the same "average" (w) ( [see 
Eq. (|SJ| below). Interestingly, as we shall see, the ISW- 
matter cross correlation breaks this degeneracy, offering 
the possibility of telling a true doomsday model from a 
model with a constant w = (w) . 

The methods presented here can obviously be used in 
a more general context, provided that there is signifi- 
cant evolution of w at low redshifts. A rather common 
assumption in phcnomenological studies of dark energy 
is to consider the simplest case of a constant w. This 
is partially motivated by degeneracies such as the one 
we just discussed above, which also occur in the angular 
spectrum of CMB anisotropics as well as in the linear 
matter power spectrum. The analysis of CMB/LSS cross 
correlations with a constant w (including w; = — 1) was 
considered in Refs. IM El IH 0. Here we shall drop 
this assumption, since the variation of w with redshift 
may provide a very exciting clue to the nature of dark 
energy, as discussed above. 

In our calculation we shall adopt a top-down approach, 
starting from the primordial spectrum of fluctuations. 
This differs from previous studies where the starting 
point is the present matter power spectrum (which is 
evolved backward in order to find its correlation with 
CMB). Our approach unifies the treatment of CMB and 
matter power spectra, and is more convenient for taking 
full account of fluctuations in the dark energy. Moreover, 
since all perturbations are evolved numerically with the 
CMBFAST code [l^, we do not resort to the frequently 
used approximate analytical expressions for the growth 
function, or the also commonly used small angle approx- 
imation [23 . The details of our calculation are reported 



in the Appendix. 

The paper is organized as follows. In Section ^ we 
describe the evolution of the universe according to the 
doomsday model. We also show that the correspond- 
ing CMB and present matter power spectra are virtually 
identical to those which are obtained in a model with a 
constant w, equal to the average (w) for the doomsday 
model. In Section llTTl we study the matter/CMB temper- 
ature cross-correlation, and show that it can be used to 
break this degeneracy. Our conclusions are summarized 
in Section Hvl 

II. EVOLUTION AND POWER SPECTRA IN 
THE DOOMSDAY SCENARIO 

We are interested in the late time evolution of our ob- 
servable universe, and so we assume a background model 
which is homogeneous and isotropic. In addition to the 
scalar field </> with linear potential the universe con- 
tains the usual radiation and matter. The dynamics is 
given by the Friedman-Robertson- Walker (FRW) equa- 
tion and the scalar field equation 

4> + ZH^ - s = 0. (4) 

Here, rimo is the fractional energy density in matter to- 
day, and _ffo is the present Hubble rate. 

In addition to the equations of motion, we need to spec- 
ify initial or boundary conditions. At early times (t — > 0), 
we expect that the scalar field is at rest and so = 0. 
This is because, in the back of our minds, we imagine 
an inflationary phase which redshifts the gradients and 
velocities of the scalar field. Hence the scalar field is ef- 
fectively homogeneous and static at i = 0. The initial 
value of the scale factor is zero as in usual Friedmann- 
Robertson- Walker evolution. The initial value of the 
scalar field is a free parameter. As mentioned in the in- 
troduction, this takes different values in distant regions 
of the universe, separated by distances much larger than 
the present Hubble radius. Finally, we require that the 
present value of the total energy density be unity. This 
is a boundary condition. These conditions can be sum- 
marized as follows: 

a(0) = 0, m^^o. m^O, H\to)^Hl, (5) 

where is the present time, defined by the requirement 
that a(to) = 1- 

The cosmological evolution following from Eqs. 13 , (ffl 
and Q has been studied by several authors 120I l2ll 
lli]. The main features of the evolution are as follows: 

• The universe starts out dominated by matter and 
hence a ~ t^/^ . At the same time the scalar field is 
essentially at rest. 
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FIG. 1: Scale factor versus time for four different slopes: s = 1 
(solid), s — 2 (short dash), s = 3 (dot) and s — 4 (long dash) 

• After some time, the matter density falls below that 
of the scalar field potential energy, and the evolu- 
tion becomes scalar field dominated. Since most of 
the energy in the scalar field is potential energy, we 
have a oc exp{Ht). 

• As the field slips down the potential, the poten- 
tial energy changes sign once 4> changes sign. With 
further slipping, a time comes when the total en- 
ergy density is zero. This epoch marks the turning 
point, where cosmic expansion changes to contrac- 
tion. From this time on, H is given by the negative 
square root of the right-hand side of Eq. 

• As the universe starts contracting, the kinetic en- 
ergy of the scalar field comes to dominate. But 
there is no stopping the contracting phase, and the 
universe rapidly arrives at the big crunch. 

These features are apparent in Fig.^ where the time evo- 
lution of the universe is represented for different values 
of the dimensionless slope 

~s = s/iVSMpH^). (6) 

One feature of this evolution that is relevant for obser- 
vational cosmology is that the equation of state for the 
scalar field changes in an unconventional manner. Very 
early on, — > and hence ~ —po where pd and 
Pd denote pressure and energy density in 0. The dark 
energy equation of state parameter is defined by 

PD 'j)^/2 - S(j) 

w = — ~ -. . (7) 

PD (/)2/2 + S0 
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FIG. 2: wiz) versus z for the four different slopes considered 
in FigEl 



Since (/) ~ at early times, w ~ — 1. At later times, 
the field starts to roll down the potential and hence the 
kinetic energy starts to play a role in po and po- This 
means that w increases with time. In Fig. |21we show w 
as a function of redshift z (Fig ures analogous to Figs.Q] 
and m can also be found in |l8ll2^ ). 

If we only varied s, keeping all other cosmological pa- 
rameters fixed, we would change the angular diameter 
distance to the last scattering surface and spoil the agree- 
ment with current CMB data. Hence it is necessary to 
suitably adjust cosmological parameters for different val- 
ues of s. The cosmological parameters in models with dif- 
ferent values of s, presented in Figs.n^ndEl depend on s. 
The simplest way to preserve the shape of the CMB fluc- 
tuations spectra on small and intermediate scales (where 
the cosmic variance is small) is to explore the so-called 
geometric degeneracy [s^ • The shape of the spectra de- 
pends mainly on two scales: - the sound horizon at 
the time of recombination, and c?a - the angular diameter 
distance to the last scattering surface. Models with same 
values of SIa/^^ and Vl\jh? have the same rg- Here Q,m 
and f2b are the total matter and baryon density fractions 
today, and h = iJo/(100 km sec~^ Mpc~^). Since the 
dark energy did not play a significant role at the time 
of recombination, changing the value of s, or w, does 
not affect r^.. The main effect changes in the properties 
of the dark energy have on the CMB spectra on small 
angular scales are due to the change in cLaj which mani- 
fests itself as a shift in the positions of the peaks in the 
angular spectra. This shift can be compensated for, with- 
out altering the structure of the peaks, by adjusting the 
value of h. As a reference model we use WMAP's best 
fit power law ACDM model 'si'l with h = 0.72, spectral 
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index n = 0.99, reionization optical depth r,. = 0.166, 
flbh'^ = 0.024, QMh^ = 0.14 and amplitude A = 0.86 (as 
defined in [^]). Given a value of s, we vary the value 



Model 


h 




s = 


0.72 


-1 


3= 1 


0.69 


-0.94 


s = 2 


0.66 


-0.81 


s = 3 


0.62 


-0.66 



TABLE I: Models considered in the paper, s is the dimen- 
sionless slope of the potential, defined in Eq. (|SJ. For each 
value of s, the dimensionless Hubble parameter h is adjusted 
so that the model reproduces the CMB peak structure ob- 
served by WMAP. The table also shows the average value of 
w for different models. 

of h, while keeping r^. A, n, flMh^ and fi^/i^ fixed, and 
find one that best reproduces the CMB spectra of the 
reference model. 

As was noted in |2^, except for the very large scales, 
CMB spectra are only sensitive to the averaged value of 
the equation of state of dark energy, (w), defined as 

da QDici)w{a) 

where the subindex Is on z and a refers to the surface 
of last scattering. In Table we show the best fit values 
of h and the corresponding values of (u')[o,2,j,] for several 
values of s. Note that the uncertainty in WMAP's esti- 
mate of h was about ±0.05 [sj. The same uncertainty 
would apply to the best fit values of h corresponding to 
models with s ^ 0. 

In the doomsday model, the equation of state param- 
eter can vary significantly at recent redshifts. This vari- 
ation is not necessarily constrained by existing analysis 
of supernova data (sS, , which by and large assumed 
a constant w. In Fig. Owe plot the effective luminosity 
as a function of redshift, as defined in for the 
s = 2 and w = —0.81, and s = 3 and w = —0.66 mod- 
els. As one can see from the Figure, the doomsday model 
with a given s and the corresponding constant w model 
have almost identical predictions in the magnitude versus 
redshift curves. Since constant w models with w > —0.8 
are disfavored at the la level js^ , the value s — 3 which 
gives an average w of —0.66 is also disfavored at the la- 
level. However, we shall still include it in our subsequent 
analysis, since it not excluded at the 2a level. 

In Fig. ^ we plot the temperature (TT) and 
temperature-polarization cross-correlation (TE) spectra 
for models in TableQ] As shown in Table I, as we increase 
s, a smaller value of Hq is needed in order to fit the CMB 
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FIG. 3: Upper panel: the luminosity, m , vs redshift plots 
for the ACDM (solid line), s = 2 (dot - short dash), s = 3 
(dot), w = —0.81 (long dash), w — —0.66 (short dash) and the 
Q,M ~ i (dot - long dash) models. Lower panel: differences 
between the m'^'^ for the models in the upper panel and the 
m"^^^ for the = 1 model. Models with the same average 
value of w are practically indistinguishable. 



spectra ^. The value h — —0.62 ±0.05 which we used for 
the s = 3 model is somewhat lower than the currently 
favored observational la region, given by ft. = 0.72 ± .08 
[37| . but still marginally consistent with it. 

The equation of state parameter w{z) obviously has 
an impact on different observables, such as CMB and 
matter power spectra. In particular, it should affect 
the integrated Sachs- Wolfe (ISW) contributions to CMB 
anisotropics. At first sight, one may think that a less neg- 
ative equation of state for dark energy would result in a 
suppression of ISW, relative to the case of a cosmolog- 
ical constant. In a universe dominated by matter, with 
pressure pm — 0, there is no time dependence of the grav- 
itational potential, and the frequency of a photon is only 
redshifted by the cosmological expansion. In this case, 
there is no late ISW effect. In the doomsday scenario, 
the dark energy equation of state is closer to that of or- 
dinary matter, and one might expect that the ISW effect 
would be smaller than in the w = —1 case. However, the 
dark energy perturbations are coupled to the dark mat- 
ter perturbations, and will also contribute to the ISW 



In principle, one could vary other cosmological parameters as 
well. However, unless one changed the model considerably, e.g. 
relaxed the assumptions of adiabaticity or scale-invariance of pri- 
mordial fluctuations, it is unlikely that one could avoid making 
Ho small in models with (w) > — 1 l3ll . 
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FIG. 4: TT and TE angular spectra for s = (ACDM) (solid 
line), s — 1 (dot - short dash), s — 2 (short dash) and s — 3 
(dot) models together with the WMAP data points. The 
value of h is chosen as indicated in Table I. Note that due to 
the geometric degeneracy, all curves look very much alike, the 
effect of a higher s being undone by choosing a lower h. 



FIG. 5: The CMB temperature anisotropy angular spectrum 
due to the late ISW effect, Cg ' , for the five models in 
Fig. |l]as well as for the model with the constant w = —0.66 
(dot - long dash line). As explained in the text, there is no 
suppression of the ISW effect due to a higher value of w. 



effect. The net result is that there is no suppression of 
the ISW effect even for values of s corresponding to aver- 
age w as high as —0.66. This is illustrated in Figs. Eland 
|6l There, we plot the angular spectrum ' and 
the auto-correlation function C^^jsw temperature 
anisotropy due to the late ISW effect, defined as 



C 



TTJSW 



(9) = {ATjsw{ni)mswin2)) 

2^+1 ^TT.ISW o / o\ 
-Co COS t'). 



E 

i=0 



47r 



(9) 



Here AT/gvy ("•) is the ISW contribution to the tempera- 
ture anisotropy in the direction h on the sky, Q is the 
angle between directions h\ and n2, and the angular 
brackets denote ensemble averaging (the expression for 
AT ISW is given in the Appendix). Note that the late 
ISW contribution of most models is in fact a bit larger 
than in the ACDM case. 

It is quite clear from Fig.^that CMB spectra alone are 
not capable of differentiating between models with differ- 
ent values of s (including s = 0) because of the geometric 
degeneracy. The effect of a larger s can be undone with 
a smaller h. (Of course, if we had a stronger observa- 
tional constraint on /i, then this would result in stronger 
constraints on s). CMB spectra alone also cannot differ- 
entiate between a model with a certain value of s and a 
model with the corresponding constant w — (w)[o,2,^] (see 
Table nj. This is illustrated in Fig. El where we plot the 
predictions of the ACDM, w = const = —0.66 and s = 3 
models. Ideally, some of the degeneracy can be removed 




FIG. 6: The ISW sourced temperature anisotropy correlation 
function C^^'^^"'(6I) for the models in Fig.^l 



by considering the matter power spectrum. Assuming 
that we had a good control over the bias, the matter 
power spectrum could in principle be inferred from ob- 
servations. In Fig. |S1 we plot the linear matter power 
spectra at z = for the w = const = —0.66 and s = 3 
models as, as well as for the s = 2 and its corresponding 
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FIG. 7: Ci vs / for ACDM (solid), s = 3 (dot) and w = -0.66 
(dash) models. 
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FIG. 8: The present linear power spectrum P{k) vs k is sen- 
sitive to (w) , but a model with a constant w has virtually the 
same P{k) as a doomsday models with the same (w). Here, 
we plot P{k) for the models in Fig.Q as well as for the s = 2 
(dot - short dash) and the w — —0.81 (long dash) models. 



w = constant — —0.81 model. Note that P{k) differs 
substantially on large scales for ACDM and the models 
with a lower average value of w. However, there is still 
an impressive degeneracy between models with the same 
(w): the curve corresponding to a constant w = —0.66 
and the curve corresponding to the model with 5 = 3 



(with (w) = —0.66) are almost identical, and similarly 
for the s — 2 and the w = —0.81 models. In the next sec- 
tion we shall concern ourselves with breaking this residual 
degeneracy. 

III. ISW AND TEMPERATURE-MATTER 
DENSITY CORRELATIONS 

In this section we show how the CMB/LSS cross- 
correlation can be used to probe the time-dependence of 
the dark energy equation of state. The cross-correlation 
is defined as 

c'^'^ie) ^ {A{h,)S{h,)) = £ ^-^crPeio) , (10) 

1=2 

where A(ni) and (5(n2) are the CMB temperature 
anisotropy and the matter density contrast along direc- 
tions ni and n2 separated by the angle 6 on the sky^ . In 
the Appendix we show that the angular cross-correlation 
spectrum Cj^"^ can be written as 

were is the primordial curvature power spectrum, as 
defined, e.g. in 38], and T^^'^(k) and Me{k) are given 

by 

T/^^ = / dTye-^^") Mk[Tj - 77o])(c$*0 - ^) , (12) 
rno 

Mi = / df] je{k[r] - Tjo]) zWg{z{r]))6{k, 77) , (13) 

where the dot denotes differentiation with respect to con- 
formal time 77, 770 is the time today, rjr is a time very early 
in the radiation era, T(r/) is the opaqueness, Wg{z) is 
the normalized galaxy selection function, (j)(k,r]), ip{k,r]) 
and 6{k,r]) are evolution functions which we define in 
the Appendix and which can be calculated numerically 
using CMBFAST 28], cg^i and are numerical coeffi- 
cients also defined in the Appendix, and j; (•) are spherical 
Bessel functions. 

The choice of the selection function Wg{z) depends on 
which large scale structure data set one wants to consider. 
Depending on the particular experiment, one also has to 
account for the possible bias between the distribution 
of the observed objects and that of the underlying dark 
matter. Our resuhs for C*^^0) for the ACDM model 
are consistent with those of ^ H, 111 1 when appropriate 
biases and selection functions are used. 

Eq (im differs from the analogous expressions in 0, |^ 
^l2^ .2^ 12^.27] as it uses the primordial curvature power 



^ The monopole and the dipole contributions depend on the choice 
of the reference frame and are not included. 
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FIG. 9: The cross-correlation breaks the remaining degenera- 
cies. In each panel, we plot a doomsday model together with 
its corresponding constant w model (these two have degener- 
ate CMB spectra and P{k)). For comparison, each plot also 
contains the fiducial ACDM model. Top two panels: C^'"^ {9) 
for the s = 2 (dot - short dash) and the w = —0.81 (long 
dash) models, as well as for the ACDM model (solid), for two 
values of Zw Bottom two plots: the same as in the top two 
plots but for the s = 3 (dot) and the w = —0.66 (short dash) 
models. 



spectrum rather than today's matter power spectrum. 
This, as explained in the Appendix, allows us to take a 
more complete account of the dark energy perturbations. 

can be particularly useful for differentiating 
between w = const and varying w models. The cross- 
correlation is sensitive to the value of w averaged over 
the range of the window function Wg{z). If w is a rapidly 
changing function of redshift, as in the case of the dooms- 
day model, then depending on where the maximum of the 
selection function is, C^'^(9) will "see" different values of 
(w) . We have calculated the cross-correlation for several 
window functions, all taken to be Gaussians of approxi- 
mately the same width as the SDSS window functions @ , 
all with the same standard deviation = 0.07 and cen- 
tered at various values of in the interval [0.1,0.8]. In 
Fig.Elwe show the plot of C^^^{e) for the S = 2 and s = 3 
models, together with their respective w — constant 
models, at two different values of z^. Corresponding an- 
gular spectra Cf^^ are shown in Fig. ^| In addition, in 
Fig. E] we plot the values of C*^'^(0.1°) as a function of 
Zw (As is seen from Fig. |51 for angular separations less 
than about a degree, the plot is insensitive to the choice 
of 0.) From Fig. E] one can see that observations focus- 
ing on redshifts in the range Zw — [0.2, 0.4] have the best 
potential of detecting the time-dependence of w. 

An interesting question is which scales give the domi- 



FIG. 10: Cross correlation angular spectra for the same mod- 
els as in Fig. |S| 
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FIG. 11: C*^^(0.1°) vs Zn, for the models in Figs. 1^ and [TI71 
ACDM (solid), s = 2 (dot - short dash), w = -0.81 (long 
dash), s = 3 (dot) and w = —0.66 (short dash). 

nant contribution to the cross-corr elation. In Fig. US we 
plot the contribution per log(fc) to the cross-correlation 
at 6* = 0.1° for the ACDM, S = 3 and w = -0.66 models 
using Zw — 0.2. Namely, we define a quantity I{k) by 

C^'^iQ.r) = J d(lnfc) I{k) . (14) 

The expression of I{k) can be deduced from Eq. (jA19|l 
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FIG. 12: I{k) vs k for the ACDM (solid), s = 3 (dot) and the 
w = —0.66 (short dash) models for z-u, = 0.2. 



of the Appendix. As one can see from Fig. ^1 I{k) has 
a broad peak around k ~ .Ol/i Mpc~^, corresponding to 
lengthscales in the range 20 — 300/i~^ Mpc. This roughly 
coincides with the peak in the matter power spectrum, 
depicted in Fig.|Sl For larger angular scales, it is still the 
same linear scales that dominate the integrand. How- 
ever, for angles 6 > 10° there is destructive interference 
between the modes and the correlation gradually disap- 
pears. 

Let us now briefly consider the prospects for obser- 
vation, and whether these will allow us to distinguish 
amongst different models. In practice, one correlates the 
CMB anisotropics with some data set {e.g a. galaxy sur- 
vey) which is supposed to trace the underlying dark mat- 
ter distribution. Of these two, the cleanest input is the 
CMB. Although very precise, current CMB data cannot 
separate the ISW contribution from the net anisotropy, 
which reduces the signal-to-noise ratio of the cross corre- 
lation. There is hope, however, that such separation may 
be possible in the future, with measurements of CMB po- 
larization towards galaxy clusters • 

A prime source of uncertainty in the matter distribu- 
tion is the bias factor. One way to proceed (which has 
been implicitly used e.g. in 5]) is the following. From the 
observed CMB autocorrelation, and following the steps 
we have described in the previous Section, the matter 
power spectrum can be normalized. Comparing this to 
the autocorrelation function of any given matter survey, 
the bias factor can be inferred. Hence, from the cross cor- 
relation of the survey with the CMB data one can infer 
the cross correlation function of the matter distribution 
with the CMB. Note in particular that the predictions 
for the matter spectra for the w = —0.66 and s = 3 mod- 



els are practically indistinguishable (see. Fig. (SJ. Hence, 
adjusting for the bias will preserve the relative difference 
between predictions for the cross-correlation for these two 
models ^ . 

Error bars in current determinations of the cross- 
correlation function H, H, IS are still too large to dis- 
tinguish between models. However, the results in have 
been obtained by using just a small fraction of the sky, 
and the situation may improve considerably with fuller 
sky coverage. Also, the uncertainties in bias may be sub- 
stantially reduced with further observations such as weak 
lensing surveys. The differences in the cross correlation 
between the models considered in this Section can be as 
large as 20%, and they have a substantial dependence 
on redshift. Hence, it seems plausible that they may fall 
within the range of detectability in a not too distant fu- 
ture. 



IV. SUMMARY AND CONCLUSIONS 

Current efforts to understand the dark energy compo- 
nent of our universe often employ dynamical scalar fields. 
The dynamics of such scalar fields lead to variation of the 
dark energy equation of state parameter, w. We have 
considered astrophysical signatures of a varying equation 
of state, 'w{z), in the context of doomsday models. There 
are two main issues that we have considered. First, can 
we tell if the average value of w is different from — I. Sec- 
ondly, what astrophysical signatures are sensitive to the 
time variation of w? 

To address both issues, we have considered autocorre- 
lations of the CMB temperature anisotropics, the matter 
power spectrum, and cross-correlation of the anisotropy 
with matter fiuctuations. The WMAP data is well repro- 
duced by models with different average value of w by suit- 
ably adjusting the Hubble parameter, and therefore the 
CMB data cannot by itself discriminate amongst differ- 
ent models. For a given CMB, the present linear matter 
power spectrum is sensitive to the average value of w{z), 
but not to its time dependence (see Fig. (HJ. Finally, the 
temperature-matter (TM)cross-correlation is sensitive to 
both the average value of w and to the time variation 
(see Figs. 151 [TUl and [TT|l . 

The TM cross-correlation requires accurate determi- 
nation of both CMB anisotropics and the matter fluc- 
tuations. Ongoing CMB observations have been very 
successful at producing data with very small error bars 



^ It should be noted that at small angles (where the cross- 
correlation with the matter power spectrum is maximal) and low 
redshifts, there is another contribution to the cross-correlation 
coming from the Sunjaev-Zeldovich (SZ) effect. Evidence for 
this contribution has already been reported in 0,13. The SZ 
contribution depends on wavelength, and a number of future 
balloon experiments (as well as the Planck satellite) will be able 
to subtract this component out. 
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and it is likely that the situation will improve even fur- 
ther in the coming years. Accurate surveys of the matter 
fluctuations are likely to be more challenging since there 
are unknowns such as the bias factor. However, since 
the CMB anisotropics and large-scale structure originate 
from the same density fluctuations, it is quite possible 
that a combination of the data can significantly reduce 
the uncertainties. We have shown that the optimal strat- 
egy for detecting a doomsday variation of w{z) is to use 
a survey in the redshift range [0.2,0.4]. 

Observations of supernovae over the next ten years are 
also likely to provide information on w{z), as discussed 
in Ref. [23| . The advantage of planned supernovae obser- 
vations are that they provide direct information on the 
Hubble expansion rate which is closely related to w{z). 
Furthermore, the data will be of good quality, with small 
error bars. However, the observations may not be too 
sensitive to rapid variation of w(z) at recent redshifts, 
the hallmark of the doomsday scenario. The TM cross- 
correlations, on the other hand, are sensitive to this time 
variation and may perhaps be used together with the su- 
pernovae observations to further constrain w{z). 
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APPENDIX A 

Here we explain details of our calculation of the cross- 
correlation. Let us define 



and 



A(n) 



Sin) 



r(n) - T 
T 



P(") - P 
P 



(Al) 



(A2) 



where r(fi) is the CMB temperature measured along the 
direction n, p{n) is the mass density along n ^, and T 



* In reality one divides the sky into pixels, with a direction n as- 
signed to each pixel, and counts the number of galaxies, N(ii) , 
inside each pixel. Then one can find the galaxy number overden- 
sity inside each pixel: (Af(n) — N)/N, which would be related to 
p(n) up to a bias factor. 



and p are the averaged CMB temperature and the matter 
density. The temperature anisotropy due to the ISW 
effect is an integral over the conformal time: 



A(n) 



dr, e--(")($-*)[(77o-77)n,r;], (A3) 



where rjr is some initial time deep in the radiation era, 
770 is the time today, $ and 5* are the Newtonian gauge 
gravitational potentials"', T{ri) is the opaqueness, which 
should, in principle, be included to account for the pos- 
sibility of late reionization, and the dot denotes differen- 
tiation with respect to 77. 

The quantity 6{n) contains contributions from astro- 
physical objects (e. g. galaxies) at different redshifts and 
can also be expressed as an integral over the conformal 
time: 

r° dz 

S{h) = y dr] — Wgizir])) (5((??o - ??)n, ??), (A4) 

where Wg{z) is a normalized galaxy selection function. 

We are interested in calculating the cross-correlation 
function 

C''^i0) = C''^i\h, ~ n2|) = (A(fii)<5(n2)) , (A5) 

where the angular brackets denote ensemble averaging 
and 9 is the angle between directions ni and £12. Let us 
introduce 



r = r]o-V 



(A6) 



The Fourier decomposition for <I'(rn,77) can be written 
as 



$(rn, 77) 



$(k,r;)e' 



(A7) 



(27r)3 

and similarly for \I'(rn, 77) and Sg(rh,r]). We can write: 



rno rva 

C^'^i9) = / dm dr,2 z{m) e-<'^'^Wg{z{m)) 



gik niT-i gik'-n2r2 



(27r)3 J (27r)3 
X ([$(k,,7i)-vi/(k,,7i)]<5(k',,72)) . (A8) 

Since the time-evolution of each Fourier mode only de- 
pends on the magnitude k = |k|, we can separate the 
directional and time dependence as: 



$(k,77) = $(k, 77^)0(^,7?) 
^'(k,r/) = *(k,?7^)7/>(fc,r7) 
5(k,77) EE (5(k,77r)5(/c,77) . 



(A9) 



Throughout this Appendix w e wo rk in the Newtonian gauge us- 
ing conventions of, e. g., Ref. |4C| . 
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Hence, we can write 3> and ^ as 



$(k, ?7) 
*(k,??) 



$(k,77^)(/.(fc,r7) 



(AlO) 



Consequently, the quantity ([3'(k, 771) — ^'(k, r]i)]6{'k', 772)) 
can be separated into the initial power spectra, which 
contain all the information relevant to the ensemble av- 
eraging, and the time-evolving part which is the same for 
all members of the ensemble: 



([$(k,77i)-vi/(k,,7i)]J(k',r72)) 
(*(k,r;,)(5(k',?7,))^(fc,r?i)J(fc',r72) 



(All) 



We take rjr to be a sufhciently early time in the radiation 
era when all modes under consideration were superhori- 
zon. Then, for adiabatic initial conditions, the growing 
mode solutions for 5, $ and "if are related to each other 
via Eil: 



C(5* 



S 



(A12) 



where Ri, = Pu/ip-y + Pi^) and pi, is the energy density 
in relativistic neutrinos. For flavors of relativistic 
neutrinos (we take — 3), after electron-positron pair 
annihilation, p^/p^ = (77V^./8)(4/ll)'*/3. This allows us 
to write: 

([$(k,,7i)-vi/(k,,7i)]<5(k',,/2)> = 
= C5*(*(k,7y^)«'(k',77^)) 

X c<E,*(/)(fc,77i) - ?/>(fc,77i) (5(fc',?72) (A13) 

From homogeneity of space it follows that 

(M'(k,r^)*(k',r;,)) = {2TT)^6^^\k + k')R^{k) , (A14) 

where Pii<(fc) is the primordial gravitational power spec- 
trum related to the more frequently used curvature power 
spectrum P-jz = 2TT'^A'^/k^ via |3a.i3aj 



invalid condition, working with the primordial spectrum 
allows us to use exact relations (|A12|) . valid deep in 
the radiation era, when all relevant modes are outside 
the horizon and dark energy fluctuations are negligible, 
and avoid the need for additional assumptions. It turns 
out, however, that, while on scales k < Q.001hMpc~^ 
and larger the fluctuations in the dark energy can be as 
large as 10%, on scales k ^ .01hMpc~^ , where the cross- 
correlation is important, the contribution of dark energy 
perturbations is rather small, less than 1%. 

Using Eqns. I|A13|) and (|A14|) we can now re- write 
Eq. lHHI) as 



q rvo rm 
^''^i(^)-25j, '^'^'J, i(^2) e--(''^)M/,(z(r72)) 



A^(fc)e'' 



k-(niri— n2r'2) 



F{k,'qi,i]2) 



(A17) 



where we have defined 



F{k,r]i,ri2) = csm c$*0(/c, ryi) - ?/'(/c, 771) S{k,r]2) . 

(A18) 

Decomposing the exponents in Eq. l|A17p into spherical 
functions and some manipulations lead to: 

C"^iO) = dr^i dm i(r/2) e^^t"^' W-g(z(7?2)) 

dk .n ,, , sm(kR) ^ , 
X / —Alik) — ^ F{k,rji,r,2) , 

(A19) 



where R = \Jr\^ r\ — 2r\r2COsQ. In addition, if one 
were to use the expression (|A19|I . one would have to sub- 
tract the monopole and dipole contributions to C*^^(0). 
This can be achieved by taking 



sin(fci?) sin(A:i?) sinfcri sinfcr2 



kR kR 
3 / sinfcri 
k'^r\r2 V kr\ 



kr\ kr2 



(A20) 



, / sinfcr2 , 

cosfcri — cosfcr2 

V 
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25^ 



(A15) 



To the best of our knowledge, in all previous litera- 
ture that contained calculations of the cross-correlation, 
it was the matter power spectrum at recent redshifts that 
was used, rather than the primordial spectrum. One can 
do that if fluctuations in the dark energy bp^ are much 
smaller than those in cold dark matter, Spcdm'- for 2: — > 0, 
one has 



$ _ ^ 



' fc2 



dm 



Spo) 



(A16) 



and usually one proceeds by assuming that 6po Spcdm 
in the equation above. While this is not necessarily an 



In practice, one wants to avoid evaluating double time 
integrals in Eq. (jA19p . A common way to reduce them 
to a single time integration is to use the so-called small 
angle {9 <^ 1) and small separation (|ri — <C ri) 
approximations |29| . These approximations were used in. 
e.g. Refs. 1^1^. One can change the integration variables 
to a; = ri — r2 and r = (ri -I- r2)/2 (or, equivalently, to 
1] ~ ivi + '72)/2) and write, in this approximation: 
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25 



dr] ie-^(") Wg{z{T])) 



dk 

- A^(fc) F{k,rj,ij) 



2'' , sin(fci?) 

dx 

-2r 



kR 



(A21) 
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where R « y/x'^ + . The integral over x can be 
evaluated analytically if one was allowed to replace the 
[— 2r, 2r] limits by [— cx),cx)]. One can assume that 2r is 
sufhciently large (r > l/kO) on relevant scales for that 
replacement to be appropriate and use 

/■°° , sin(fci?) TT ^ 
J kR " k^^^''^^ (^22) 

to obtain the following form: 

X A^ik) Jo{ke[m-v]) F{k, f],fi) . 

(A23) 

While the approximate expression (|A23|I can be useful 
for analytical estimates, it gives an error of order 2 — 4% 
on the scales of interest and we shall not resort to it. 

Instead of evaluating the expression l|A19p directly, 
from computational point of view, it is advantageous to 
decompose it into Legendre series, compute the individ- 
ual coefficients of the decomposition, and then sum the 



series. Namely, Eq. (|A19|) can be written as 

C''^{9) - ^IT^^f "^^^(cos 0) , (A24) 

1=2 

where we do not include the monopole and dipole terms 
in the sum, and where can be written as 

Cr = 4vrA I ^ T/^H'(fc) M,{k) , (A25) 

with functions Tg^^{k) and Mi{k) defined as 
rVo 

r/^'^= / e~^Wdryj,(fc[77-77o])(c3.*,^-^)(A26) 
Me = cs^ / drj j,(k[f^ - ,^^])zWg(7M)hk, v) ,(A27) 

where are spherical Bessel functions. One can use 
CMBFAST 28], with minor modifications, to compute 
functions Tf^'^ and and to normalize A^(fc). 
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